_ ^STER  COPY 


AD- A247  528 


KEEP  THIS  COPY  FOR  REPRODUCTIOM  PtjRPOSES 


Form  Approvtd 
0MB  No.  0704-0188 


(VTION  PAGE 


•  1  hour  o«r  roipofltt,  mclud(n9  ih«  limt  lof  rtvitwing  inttrurtionft,  >«sr(hmg  tKiltmg  Oita 

ng  tha  (ollactton  of  informaiioo.  Sand  <ommaAt»  regarding  thit  burdtrt  aitlmatt  or  arty  other  aipect  of  thii 
^  .  .  to  «A/a%hington  Headduartert  Servicei,  Directorate  for  information  Ootrationi  and  Mportt,  1215  Jtfftrvon 

uavit  Highway.  Suite  UQ4.  Arlington.  VA  22202<4302.  and  to  me  Office  of  Management  and  Budget.  Paperwork  Reduction  Project  (0704*0  IBB),  Wathington.  DC  20503< 


1.  AGENCY  USE  ONLY  (Lfive  blank) 


2.  REPORT  DATE 

January  16,  1992 


3.  REPORT  TYPE  AND  DATES  COVERED 


4.  TITLE  AND  SUBTITLE  ^lonte  Carlo  Simulations  of  Random  Rough 
Surface  Scattering  with  Finite  Element  and  Finite 
Difference  Methods 


6.  AUTHOR(S) 


Leung  Tsang  And  Chi  Hou  Chan 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


9.  SPONSORING /MONITORING  AGENCY  NAME(S)  AND  ADDRESS(tS) 


U.  S.  Army  Research  Office 
P.  0.  Box  12211 

Research  Triangle  Park,  NC  27709-2211 


10.  SPONSORING /MONITORING 
AGENCY  REPORT  NUMBER 


11.  SUPPLEMENTARY  NOTES 

The  view,  opinions  and/or  findings  contained  in  this  report  are  those  of  the 
author (s)  and  should  not  be  construed  as  an  official  Department  of  the  Army 
position,  policy,  or  decision,  unless  so  designated  by  other  documentation. 


12a.  DISTRIBUTION /AVAILABILITY  STATEMENT  12b.  DISTRIBUTION  CODE 

Approved  for  public  release;  distribution  unlimited. 


13.  ABSTRACT  (Maximum  200  words) 

Monte  Carlo  simulations  of  scattering  of  waves  by  one-dimensional  random 
rough  gold  surfaces  at  optical  frequencies  are  investigated  by  the  use  of  the 
finite-element  method.  Both  TE  and  TM  plane  wave  incidences  are  considered. 
Backscattering  enhancements  are  observed  for  all  the  cases  studied  in  this 
report . 


DEFENSE  TECHNICAL  INFORMRTION  CENTER 


9a05939 


14.  SUBJECT  TERMS 

Random  rough  surface,  Monte  Carlo  simulations,  backscattering 
enhancement,  finite-element  method 


15.  NUMBER  OF  PAGES 

4 


16.  PRICE  CODE 


17.  security  CLASSIFICATION  118.  SECURITY  CLASSIFICATION  119.  SECURITY  CLASSIFICATION  I  20.  LIMITATION  OF  ABSTRACT 


OF  report 

UNaASSIFIED 


NSN  7540-01 -.’BO-SSOO 


OF  THIS  PAGE 

UNCLASSIFIED 


OF  ABSTRACT 
UNCLASSIFIED 


Standard  Form  298  (Rev  2-89) 

PtrunM  by  Mb  JJS-'a 


FINAL  REPORT 


January  16,  1992 


1 .  ARO  Proposal  Number:  29276-GS 


RTV;  tjHkkl 
VTic  rmi 
iJfWji.'ino"jiu  art 
Josclf  loatlo*! 


2.  Period  Covered  by  the  Report:  July  15, 1991  to  November  14, 1991 


3.  Title  of  Proposal:  "Monte  Carlo  Simulations  of  Random  Rough  Surface  Scattering  with 

Finite  Element  and  Finite  Difference  Methods"  i-ttbU-itr  Co*«fG 


4.  Contract/Grant  Number:  DAAL03-91-G-0217 

5 .  Name  of  Institution:  University  of  Washington 

6.  Authors  of  Report:  Leung  Tsang  and  Chi  Hou  Chan 


r:* 


(Afilil  HMfi/ar 


lat  1 

1 

Kpflfllal 

1.  Scientific  Personnel  Supported  By  this  Project: 
Leung  Tsang,  Principal  Investigator 
Chi  Hou  Chan,  Co-Principal  Investigator 
Charles  Mandt,  Research  Assistant 


Summary  of  the  Research  Findings 

The  earlier  formulation  of  Monte  Carlo  simulations  of  scattering  of  waves  by 
perfectly  conducting  random  rough  surfaces  based  on  the  finite-element  method  and 
periodic  boundary  condition  is  extended  to  investigate  scattering  from  one-dimensional 
random  rough  gold  surfaces  at  optical  frequencies.  In  this  report,  two  cases  are 
investigated.  For  Case  1,  the  random  surfaces  have  a  correlation  length  of  3.099X  and  a 
root-mean-square  height  of  1.6927X.,  and  the  refractive  index  of  0.312  +  i7.93  is  used. 
For  Case  2,  these  numbers  are  1.0525X,  0.57 49A.,  1.958  +  i20.7,  respectively.  Both  TE 
and  TM  plane-wave  incidences  at  0®,  10®,  and  30®  are  considered.  The  finite  conducting 
surfaces  are  modeled  as  lossy  dielectric  rough  surfaces.  To  reduce  the  computational 
domain,  for  the  TE  case,  we  insert  a  perfect  electric  conducting  (PEC)  wall  which  has  the 
same  rough  surface  profile  at  a  few  skin  depths  below  the  air-dielectric  interface. 
Similarly,  for  the  TM  case,  a  perfect  magnetic  conducting  (PMC)  wall  is  used. 

In  our  formulation,  the  scattered  field  slightly  above  a  reference  plane  at  a  small 
distance  (0.5>.)  above  the  maximum  height  of  the  rough  surface  is  expanded  in  terms  of 
Floquet  modes.  The  region  bounded  by  this  reference  plane  at  the  top,  the  conducting 
plane  at  the  bottom,  and  the  periodic  boundary  condition  on  the  two  sides  is  calculated 
using  the  finite  elements  with  the  field  at  the  reference  plane  being  specified  as  cosine  or 
sine  functions.  We  call  these  solutions  finite-element  modes.  There  is  no  one-to-one 
correspondence  between  the  FEM  modes  and  the  Floquet  modes.  In  fact,  each  of  the  FEM 
modal  solutions  in  the  region  very  close  to  the  rough  surface  incorporates  all  the 
propagating  modes  and  a  large  number  of  evanescent  modes.  Most  of  these  evanescent 
modes,  however,  decay  to  zero  at  the  reference  plane.  For  each  of  the  finite-element 
modes,  we  solve  a  Helmholtz  equation  in  which  we  use  a  discretization  of  6  nodes  per 
linear  wavelength  for  all  the  calculations  reported  here.  The  FEM  calculation  is  very 
efficient  with  the  use  of  a  sparse  matrix  solver.  The  amplitudes  of  the  Floquet  modes  are 
calculated  by  matching  the  field  and  its  normal  derivative  represented  by  the  Roquet  modes 
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and  linear  combinations  of  FEM  modes  on  the  two  sides  of  the  reference  plane.  The  fast 
Fourier  transform  (FFT)  algorithm  has  been  used  in  the  process  to  improve  numerical 
efficiency.  Extensive  numerical  experiments  show  that  only  four  to  six  evanescent  modes 
are  needed.  The  maximum  error  in  the  power  conservation  is  less  than  1  %  when  the  lossy 
medium  is  replaced  by  a  lossless  one.  There  are  two  propagating  Floquet  modes  for  each 
wavelength  of  the  surface  length.  Therefore,  a  full  matrix  equation  of  the  order  of  2 
[integer  (L  /  \)]  +  Ne  +  I  needs  to  be  solved,  where  Ne  is  the  number  of  evanescent  modes 
and  L  is  the  surface  length.  In  the  calculation  of  surface  length  of  SOX,  we  have  used  101 
propagating  modes  and  4  evanescent  modes,  a  total  of  105  Floquet  modes.  The  same 
number  of  FEM  modes  are  used.  In  contrast,  for  the  spectral  domain  extended  boundary 
method  in  which  the  field  is  expanded  right  on  the  rough  surface,  as  many  evanescent 
modes  as  propagating  ones  are  often  required  for  the  large  surface  height  considered  in  this 
study.  On  the  other  hand,  a  full  matrix  of  the  order  1000  x  1000  is  required  for  the  integral 
equation  method  for  this  two-medium  problem  when  a  discretization  of  10  per  linear 
wavelength  is  used  for  the  rough  surfaces  under  consideration. 

In  the  numerical  results  to  follow  we  present  the  radar  cross  sections  of  the  rough 
surface  for  both  TE  and  TM  incidences  at  0®,  10®,  and  30®  at  the  wavelengths  of  1.152  trm 
(Case  1)  and  3.392  pm  (Case  2).  To  investigate  the  effect  of  the  surface  length  L  to  the 
radar  cross  section,  two  surface  lengths,  namely,  L  =  30X  and  50X  are  us  'd  for  Case  2  with 
the  TM  incidence  at  30®.  Figure  1  shows  that  the  RCS  exhibits  two  peaks,  one  at  the 
backscattering  direction  and  the  other  at  the  specular  direction,  when  the  surface  length  is 
small.  In  Figure  2,  however,  the  specular  peak  is  suppressed  while  the  backscattering  is 
enhanced  when  we  increase  the  surface  length  to  SOX.  For  this  two-medium  problem,  due 
to  limited  computer  resources,  we  chose  all  our  surface  lengths  to  be  SOX.  Figures  3  to  14 
show  the  RCSs  for  all  the  12  cases.  We  summarize  our  findings  as  follows:  (1)  both  the 
TE  and  TM  incidences  show  backscattering  enhancement  at  all  three  incident  angles;  (2)  it 
is  also  found  that  while  the  TE  and  TM  responses  are  vei7  similar  in  Case  1,  TE  incidence 
yields  stronger  backscattering  enhancement  in  Case  2,  in  contrast  to  the  scattering  from 
surfaces  with  small  roughnesses  where  backscattering  enhancement  is  observed  only  for 
the  TM  waves;  (3)  of  all  the  numerical  simulations,  the  lossy  surfaces  absorb  more  power 
for  the  TM  incidences  than  that  of  the  TE;  and  (4)  the  absorption  rate  decreases  with  the 
incident  angles. 


Case  2 

200  Kcall’talluns 
I,  s  30,2  \Vavclcii||lhs 


Case  2 

200  Kealizalloiis 
L  a  S0.2  Wavelengths 


SCATTERING  ANCLE  (degree) 
rigurc  2.  Nonnalized  cross  section  versus  Ksliering  angle. 
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Figure  3.  Nonnalizcd  crou  section  versus  scsitering  snglc. 
TE  incidence  si  0  degree. 
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Figure  S,  Nomislized  cross  section  versus  Ksllering  snglc. 
TE  incidence  at  30  degree. 


Case  2 

200  Keallzalions 
1.  >30,2  Wavelengths 


SCATTEKINC  ANCLE  (degree) 

Figure  7.  Nonnslized  erots  section  versus  Ksltering  snglc. 
TB  incidence  at  10  degree. 
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Figure  4.  Nomtsliied  cross  section  veisus  scattering  angle. 
TE  incidence  si  10  degree. 
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Figure  6.  Noramlized  cross  section  versus  scattering  angle. 
TE  incidence  at  0  degree. 
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111  incidence  at  30  degree. 
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TM  Incidence  alio  degree,  1 M  Incidence  al  30  degree. 


